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ARRANGEMENT FOR FORECASTING TIMELY 
COMPLETION OF A TASK BY A RESOURCE 

Technical Field 

This invention relates to work management systems in general 
5 and to communications-serving and -distributing systems in particular. 

Background of the Invention 

Work management systems, including call management, 
workflow automation, work distributor, and automatic call distribution 
(ACD) systems, distribute work items— whether tangible or intangible, and 

10 referred to herein as tasks— for servicing to resources, such as processing 
machines or call-center agents. In order to operate efficiently, work 
management systems must anticipate when resources will become 
available to serve a new task, so that they can have a task ready for a 
resource as soon as the resource becomes available and preferably no 

15 earlier. Otherwise, tasks such as calls must wait in queue for resources to 
become available, to the dissatisfaction of the callers, or resources sit idle 
while awaiting tasks to become available, to the detriment of operating 
efficiency of the system. 

Work arrangement systems often do not use an accurate 

20 forecast of the future availability of presently unavailable resources. In 
these kinds of systems, resources are presently unavailable because they 
are performing tasks that are designated to be uninterruptible. That is, a 
resource can serve only one task at a time and cannot start a new task 
until the present task is completed. Additionally, only one resource can 

25 usually be assigned to a particular task. Illustrative such tasks involve live 
clients in transactions, such as telephone calls. Background tasks may be 
dynamically interspersed amongst them. Resources that are presently 
handling tasks are rendered presently unavailable. However, each 
unavailable resource can be expected to complete its task within a future 

30 time interval with a determinable probability. This expectation for the 
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"arrival" of a resource can be based upon its time-in-state and the 
particular operational characteristics of the type of task, which the 
resource is serving. 

Typically, automatic systems forecast the arrival of new tasks 
5 better than they forecast the arrival of resources to serve them. For 
example, one known call-management system does not use dynamic 
forecasts at all. Instead, its prediction algorithm triggers the initiation of 
new outgoing calls (tasks) for each agent (resource) after the agent has 
been in a 'work', 'record update", or 'wrap-up' state for a relatively fixed 

10 amount of time. The system does not compute a probability of task 

completion for agents, either collectively or individually. This results in low 
agent utilization, and in high client nuisance rates due to clients' time 
spent waiting in queue and client-abandoned calls. 

A task-acquisition subsystem should strive to keep resources 

15 busy serving tasks at all times. For example, an outbound predictive 

dialer should initiate outbound calls based upon the aggregate availability 
of agents. The time horizon for this determination should be close to the 
amount of time that it takes to have an incremental outbound call attempt 
answered by a human being (i.e., to obtain a live-party answer). Then 

20 outbound calls can be initiated to generate a demand for service that 
appropriately reduces the gap between the supply of agents and the 
demand for them. 

An automatic call distributor or work distributor dispatches a 
resource to serve lower-priority tasks when the supply of resources 

25 exceeds by some selected amount the demand for service by higher- 
priority tasks. In this case, the relevant time horizon for a forecast is 
typically close to the amount of time that the resources could be assigned 
to lower-priority tasks before becoming available again to service other 
tasks. 
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Summary of the Invention 

This invention is directed to solving these and other problems 
and meeting these and other requirements of the prior art. The invention 
provides an effective measure of an individual resource's contribution to 
5 the supply of available resources. It determines the contribution that a 
resource makes to resource availability. The aggregate resource supply 
may then be taken as the sum of the availabilities of the individual 
resources. 

According to one aspect of the invention, a work-management 

10 method determines a probability of availability at a future point in time of a 

plurality of resources, combines (illustratively sums) the probabilities of the 

individual resources to obtain a number, and schedules for the future point 

in time no more than the number of new tasks for servicing by the plurality 

of resources. According to another aspect of the invention, a work- 

15 management method determines an amount of time f that a resource has 

been servicing a task by now, determines a probability F(t+h) of the 

resource servicing the task to completion within a total amount of time t+h, 

determines a probability F(t) of the resource completing servicing the task 

by now (within time Q, determines a probability P that the resource will 

20 complete servicing the task within an amount of time h from now as 

F(t + h)-F(t) . . . , 

— , and in response to P schedules another task for 

l-F(t) 

servicing. l-F(t) is illustratively the equivalent of a determined probability 
of the resource not completing servicing the task by now (within time t). 
Preferably, the determining steps are performed for a plurality of 
25 resources, a number of the resources that will likely have completed 
servicing their respective tasks within the amount of time h from now is 
determined by summing the probabilities P that have been determined for 
the individual ones of the plurality of resources, and the number is used to 
schedule new tasks for servicing. Preferably, no more than the number of 
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new tasks are scheduled for the amount of time h from now for servicing 
by the plurality of the resources. 

This method is significant because it enables the realization of 
high service levels and high utilization of resources simultaneously. It can 
5 be used to improve customer service and reduce the burden on 

supervision in work-processing centers, including call centers and contact 
centers. A predictive dialer may use this method to avoid annoying call 
recipients with outbound calls for which no agent will be available to 
service them. A workflow system can use this method to ensure that it 

10 does not assign an available resource to low-priority work when the 
resource should be dedicated to higher-priority work. 

While the invention has been characterized in terms of 
method, it also encompasses apparatus that performs the method. The 
apparatus preferably includes an effecter— any entity that effects the 

15 corresponding step, unlike a means— for each step. The invention further 
encompasses any computer-readable medium containing instructions 
which, when executed in a computer, cause the computer to perform the 
method steps. 

These and other features and advantages of the invention will 
20 become more apparent from the following description of an illustrative 
embodiment of the invention considered together with the drawing. 

Brief Description of the Drawing 

FIG. 1 is a block diagram of a call center that includes an 
illustrative embodiment of the invention; 
25 FIG. 2 is a functional diagram of a resource task-completion 

forecaster of the call center of FIG. 1 ; and 

FIG. 3 is a functional diagram of a task scheduler of the call 
center of FIG. 1. 
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Detailed Description Of An Illustrative Embodiment 

Reference will now be made in detail to the illustrative 
embodiment of the invention, which is illustrated in the accompanying 
drawing. While the invention will be described in conjunction with the 
5 illustrative embodiment, it will be understood that it is not intended to limit 
the invention to this embodiment. On the contrary, the invention is 
intended to cover alternatives, modifications, and equivalents, which may 
be included within the invention as defined by the appended claims. 
FIG. 1 shows an illustrative call center 100 for servicing 

10 inbound and outbound calls. The word "call" is used herein generically to 
mean any communication or request for expeditious service. Call 
center 100 comprises an automatic call distributor (ACD) 104 that 
interconnects agent positions 130-140 via calls with the outside world to 
which it is connected by communications trunks 102. ACD 104 includes a 

15 switching fabric 1 16 that selectively interconnects trunks 1 02 with 
communications lines 106 that extend to agent positions 130-140. 
ACD 104 is a stored program-controlled apparatus that operates under 
control of a processor 112 that obtains and stores data in, and executes 
stored programs out of, memory 1 10 or any other computer-readable 

20 medium. Data in memory 110 include historical and operational data of 
ACD 104 and agents 156, which are stored in a call management system 
(CMS) 120 database. Processor 112 controls operation of switching 
fabric 1 1 6 and of a dialer 1 14 that generates outgoing calls on trunks 1 02 
through switching fabric 116. Each agent position 130-140 includes a 

25 terminal 152, such as a personal computer, and a voice communications 
device, such as a telephone or a headset 154, for use by an agent 156. 
As described so far, call center 100 is conventional. 

For purposes of the following discussion, a call, whether 
incoming or outgoing, constitutes a task to be served, and an agent 

30 position 1 30-1 40 that is presently staffed by an agent 1 56 constitutes a 
resource for serving tasks. According to the invention, memory 1 10 of 
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ACD 110 includes a resource task-completion forecaster 122 program and 
a resource scheduler 124 program. Forecaster 122 forecasts a probability 
of completion of servicing a call by a certain time by an agent who is 
presently assigned to service the call. The forecast for each agent 
5 indicates a probability that the agent will complete servicing the presently 
assigned call, and consequently will become available to service another 
call, by a specified future point in time. Forecaster 122 thus provides a 
measure of contribution of individual agents 156 to the aggregate supply 
of available agents for use by resource scheduler 124 in automatically 
10 pacing the generation of outbound calls by dialer 1 14 and in automatically 
allocating agents 156 to calls. The configuration and operation of 
forecaster 122 is shown in FIG. 2 and is described in conjunction 
therewith. 

Forecaster 122 performs the sequence of nine basic steps 

15 shown in FIG. 2. It is illustratively executed for each agent 156 at any time 
that ACD 104 could have a need for an agent 156 or a need for work at 
some future point in time (a horizon interval h). Forecaster 122 could 
execute periodically and whenever a call or an agent changes state. 
Upon its invocation, at step 200, forecaster 122 initially determines the 

20 type of call that the agent 1 56 whose task-completion forecast it is 

attempting to compute is presently serving, at step 202, and from the call's 
record determines the amount of time that the agent 156 has already 
spent servicing this call, at step 204. The time determined at step 204 is 
called "time-in-state" t. 

25 Forecaster 1 22 then retrieves from the CMS 1 20 the task- 

completion characteristics for this type of call, at step 206, and fits the 
task-completion statistics into a lifetime closed-form cumulative-probability 
distribution F, at step 208. Forecaster 122 then numerically evaluates the 
cumulative-probability distribution F for the time-in-state fas F(t), at 

30 step 210, and also numerically evaluates the cumulative probability 
distribution for a time comprising the sum of the time-in-state t and the 
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selected time horizon h as F(t + h), at step 212. These steps require prior 
execution of steps 204 and 208. Forecaster 122 then computes the a- 
priori probability of completion of servicing the call within the horizon 
interval as R=F(t + h)-F(t), at step 214, and also computes the a-priori 
5 probability of non-completion of servicing the call before the horizon 
interval expires as Q = l-F(t), at step 216. These steps require prior 
execution of steps 210 and 212. Finally, forecaster 122 computes the 
conditional probability of completion of servicing the call within the forecast 

horizon as P = — at step 218. Step 218 requires prior execution of 

10 steps 214 and 21 6. Forecaster 122 then ends its operation, at step 220. 

Step 202 uses the attributes of the call to classify it. When an 
agent 156 is servicing a call, typically the call is already classified to a 
great extent by conventional means. Such classification is typically 
needed by ACD 104 for the purpose of routing the call and distributing it to 

15 agent positions 130-140. The call can be differentiated by the channel, 
the relationship, the customer type, the destination directory number, etc. 
Examples of the type of channel are inbound voice call, voice message, 
outbound telephone call, e-mail, FAX, letter, and World Wide Web hit. 
The relationship could be existing customer, prospective customer, or 

20 unknown. The existing-customer type could be preferred, normal, or 

marginal. The destination directory number can determine the motivation 
for the call or the expected type of transaction, such as inquiry, loan 
request, order, sales, post-sales customer service, etc. Additionally, the 
customer or prospect record and the skill attributes of the agent can 

25 provide a myriad of attributes that distinguish the calls. Generally, the 
more fully that a call is classified, the less variance is there in the work 
statistics for completion of the servicing of the call. 

Step 204 generally involves retrieving a clock time that was 
stored in the call record when the call was distributed to the agent, and 

30 then subtracting that clock time from the present clock time. This time-in- 
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state is important because the outlook for completion of work generally 
changes as the work progresses. In contrast, most commonly-used 
queuing models assume that the outlook for completion of a transaction is 
invariant with respect to the amount of time already spent working on it. 
5 This step does not depend upon any previous steps. 

In this illustrative embodiment, the task completion statistics 
that are retrieved at step 206 comprise at least the mean of the work time 
(i.e., first moment of the work time) and the variance of the work time (i.e., 
second moment of the work time) that either this particular agent or all 
10 agents on average have historically spent on servicing this type of call. 
The first moment (mean) and the second moment (variance) are 
conventional statistical measures. The sample mean is typically defined 
as the sum of the values of the samples divided by the number of 
samples: 

1 " 

15 M = — ! xi 

i = l 

while the sample variance is typically defined as the sum of the squares of 
the differences between the values of the samples and the mean divided 
by the number of samples less one: 

ct 2 =— ^U-M) 2 
n — l , = , 

20 Recent samples can be weighted more heavily by employing 

moving averages or exponential smoothing in the definitions of the mean 
and variance. The mean and variance could be determined by time and 
motion studies, or other theoretical models of agent or client behavior. 
Alternatively, histograms could be used to describe the distribution without 

25 benefit of descriptive statistics. Typically, these statistics would have been 
generated from samples of historical results for the subject type of call. 

Step 208 involves fitting the statistics that were retrieved at 
step 206 to a classical closed-form cumulative-probability distribution, or 
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alternatively representing the statistics in a histogram, in order to 
determine values for a dispersion parameter (a) and a parameter of 
central tendency (p). For ease of computation, the illustrative embodiment 
uses a classical two-parameter Weibull distribution, which is a 
5 generalization of the exponential distribution. It gives the distribution of 
lifetimes of objects. The Weibull distribution is defined as 

-{tlp) a 

F(f) = l-* 

which gives the probability of the event occurring before time t The mean 
is defined as // = /? T(l + a' ) . 

10 The variance is defined as a = fi 1 [r(l + 2a' 1 ) - r 2 (l + a' 1 )]. 
The skewness is defined as 

2T 3 (1 + a~ l ) - 3r(l + a~ l )T(l + 2a~ l ) + T(l + 3a' 1 ) 

vl = — 

[r(i + 2a" 1 )-r 2 (i + a~ 1 )J /2 
and the kurtosis is defined as 

r(l + 4a 1 ) - 4r(l + 3a 1 )T(l + a' 1 ) + 3T (1 + 2a~ l ) 
v2= -6 

[r(i + 2a" 1 )-r 2 (i + cT 1 )] 
15 where t e [0 ? <>o] ?e > j s the base of natural logarithms, a is the dispersion 
parameter, (3 is the parameter of central tendency, and r is the classical 

"gamma" function defined as T(x) = | y x ~ x e' y dy , where x> o. 

The following steps describe a method of numerically 
evaluating the classical gamma function. This method is simpler, quicker, 
20 but less accurate than other known methods. 

1 . Assign the input argument's value to x. 

2. Enforce a positive argument by x <r- max(jc, 0.0001). 

3. If x < 1 , then z <- x + lto assign a value to the input parameter of 
the core method for arguments greater than 1 . 
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4. If x > 1, then z <- x to assign a value to the input parameter of the 

core method for arguments greater than 1 . 

5. y <r- z + 4.5 . 



6. w<—ln 



" 76.180091729406 86.505320327112 24.01409822223 1.23173951614' 
1 + + 



z z+1 z+2 z+3 

5 where In is the natural logarithm function, with a base of e . 

7. w<- w+(z-0.5)log(y)-y+ln( < -j2jvy This is the natural logarithm 

of the gamma function for an input argument not less than 1 . 

8. If x 56 z , then w <r- w- ln(x) to get the natural logarithm of the 

gamma function in the case where the input argument is less 
10 thanl. 

9. Compute the gamma function value from the natural logarithm of 
the gamma function value by r <- e w . 

Given the mean ju and the variance a 1 about the mean, one 
method of fitting statistics to the two-parameter Weibull distribution 
15 determines an approximation for the parameters a > 0 and fi > 0 for the 
cumulative probability Weibull distribution as follows. The method 

determines the value of a by a = — . It determines x by iteration as 



x, 



follows. First, x 0 =1 ; then, for each succeeding +1* iteration 

x,[g(2x,)-2g(x,)] 



g(2x l )-g{x i ) 



(o 2 } 

20 where g(x)=log(r(2x+l))-21og(r(x+l))-log — +1 . Normally three or four 



V" J 



1 1 

iterations are sufficient, i.e., a = — or — . Then the method 



'4 



determines from the mean and a by p-- 

r - + i 
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For example, assume that the mean // is 1 10 seconds and the variance 
cr 2 about the mean is 18150 seconds squared. Then 

log[t7 2 //z 2 +l]=log(2.5) = 0.9162907. Withx=1, T(x+l) = T(2) = l and 

T(2x+1) = r(3) = 2; and g(x) = log(2) - 2 log (1)- 0.9 162907 = -.22314. 
5 T(4x + 1) = T(5) = 24, so g(2x) = log(24) - 2 log (2) - 0.9162907 = 0.8754687. 
Consequently, the next value for x is 

x, (0.8754687 - 2 * (-0.22314)) (l)( 0 .8754687 + 0.44628) 

*i = = — 1 - = 1.2031 . 

0.8754687 - (-0.22314) 1 .098609 

Then g(x,) = 0.014585116 and g(2x,) = 1.3564625 giving the next x as x 2 =1.215913. 
g(x 2 ) = 0.001 112875 and §(2^) = 1.387193 giving the next x as X3 =1.216887. 

1 1 

10 Taking x 3 to be x, the method gives a=-= = 0.82178047. 

x 1.216887 



Finally, p = — ^ = _ ]}°_ = ^ = 98.914909 seconds , 

r 



+ T(2.216887) 1.112066933 



\ a J 

So, the method determines the fit to the Weibull distribution to be 

F(t) = 1 - e -( f/98 -914909)0-82178047 

A variety of means can be employed to perform the fit. Typically, this 
15 involves some iteration, and errors can be expected. One can reduce the 
magnitude of errors by using more parameters to describe the distribution. 
Step 208 requires the previous execution of step 206. 

At step 210, the probability distribution Fis evaluated for time 
value t. This step requires the previous execution of steps 208 and 204. 
20 At step 212, the horizon interval h is selected based on when 

in the future the system expects to have a need for an agent 156. For 
example, if it is known that it typically takes x seconds from the time that 
dialer 114 dials a call to the time when the call is answered, and the 
system is trying to determine whether or not to cause dialer 1 14 to dial a 
25 call at this instant, the value of h may be selected to be x. 
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Step 214 gives the a-priori probability of completion of 
servicing the call during the horizon interval without consideration of the 
amount of time that has already been spent servicing the call without 
completing it. This is the outlook for completing the servicing of the call in 

5 the horizon interval assuming that servicing of the call has not yet started. 
In practice, this value must usually be attenuated or subject to a floor in 
order to avoid unrealistically low values that are not supported by sample 
data from historical experience. 

Step 216 gives the a-priori probability that the agent has not 

10 completed servicing the call by the time-in-state. This value requires the 
same attenuation as is applied in step 214. 

Step 218 computes the conditional probability of completion in 
the horizon interval, given the agent's tenure on servicing the call (the 
time-in-state). When an attenuation scheme is applied to steps 214 

15 and 216, the result of step 218 should approach zero asymptotically as the 
time-in-state increases. Thus, a proper attenuation scheme effectively 
discounts the "arrival" of agents (i.e., agents becoming available) with 
"outlier" times-in-state. "Outlier" times are times that are grossly out-of- 
normal. 

20 The basic functionality of resource scheduler 124 is shown in 

FIG. 3. Upon receiving the probabilities P for each agent 156 calculated 
by forecaster 122, at step 300, it combines those probabilities, preferably 
by summing them, at step 302, to obtain a number. The sum of the 
probabilities P determined for all agents 156 represents the expected 

25 number of available agents 156 within the forecast horizon h. Resource 
scheduler 124 then uses the number to schedule new tasks (calls) for 
agents 156, at step 304. Preferably, scheduler 124 schedules no more 
than that number of new calls to become available for servicing by the 
agents for the forecast horizon. Illustratively, scheduling that number of 

30 new calls involves directing dialer 1 14 to place as many outbound calls as 
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are expected to result in that number of calls being answered by a person. 
Resource scheduler 124 then ends its operation, at step 306. 

The invention is useful for the automation of resource 
allocation and various dispatching operations, especially in an 

5 environment of differentiated tasks and differentiated resources where 
control of expeditious service is required for some of the tasks. For 
example, a predictive dialer may use this method to avoid annoying call 
recipients with outbound calls for which no agent will be available to 
service them. Workflow systems can use this method to ensure that they 

10 do not assign available resources to low-priority work when the resources 
should be dedicated to higher-priority work. The invention can thus 
facilitate the integration of back-office work with front-office call work, 
including the integration of outbound call systems with workflow systems. 

Of course, various changes and modifications to the illustrative 

15 embodiment described above will be apparent to those skilled in the art. 
For example, the tasks could have possible outcomes of failure or success 
in that the agent could become unavailable with some probability, the 
tasks could involve serial steps each of which has statistics facilitating an 
arrival forecast, the forecast could depend upon the statistics of the path 

20 of the automated script taken by the agent, the forecasting method could 
apply to the generation of work subject to lead time uncertainties, the 
scheduler could schedule work based on criteria for answer delay 
including service levels, the scheduler could schedule work based upon 
utilization of resources, a forecast could be performed for several time 

25 horizons h and the scheduler could employ necessary or sufficient criteria 
for each forecast, the scheduler could employ this method conditionally 
based on normal conditions while under abnormal conditions employing 
another method, the forecasting could employ handicaps for agents so as 
to get a better forecast for each individual agent, the statistics and 

30 forecasting could be differentiated by agent type and call type, individual 
agent statistics could be updated asynchronously and independently from 
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the generation of the aggregate arrival statistics, the frequency of updating 
an individual arrival probability could depend on the state of the work on 
the call, the update of each call could be programmed for only certain 
times in state (e.g. 0, 5, 20, 40, 60, 80, 100, and 150), or the forecast 

5 horizon h could adjust dynamically to systematically provide a "worst- 
case" forecast. Such changes and modifications can be made without 
departing from the spirit and the scope of the invention and without 
diminishing its attendant advantages. It is therefore intended that such 
changes and modifications be covered by the following claims except 

10 insofar as limited by the prior art. 
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